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ABSTRACT 


The non-linear equations for an incompressible turbulent 
atmosphere are integrated numerically to investigate the sea-—breeze 
circulation. Solutions are presented for various initial states of 
the atmosphere. In addition, liquid water and water vapour are 
incorporated to permit the study of the effect of the sea-breeze 


circulation on a coastal fog bank. 
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CHAPTER 1 


THE SEA-BREEZE 


DAL introduction 

In coastal areas, during conditions of strong insolation and 
nocturnal cooling, there is a marked local diurnal variation in the 
surface winds. During the day, a wind (called a sea-breeze) blows from 
sea to land and during the night, a reverse flow from land to sea (a 


land-breeze) is observed. 


de2 Déescrrocvrom or “the Crremlation 


During the day, along the sea coasts or shores of large lakes, 
a temperature gradient is established as a result of the differential 
warming of the land and sea by the same heat source. According to 
Defant (1951), the solar energy incident upon the water is distributed 
through a deep layer by conduction or turbulent mixing, resulting ina 
small temperature rise. The land, on the other hand, has only a shallow 
layer heated by the sun; hence, a relatively large temperature rise occurs. 
Observations by Hsu (1967) show that the surface temperature gradient is 


concentrated along the coast. 


The heated land causes the temperature of the air above it to 
rise as a result of turbulent diffusion and thermal convection. In the 
warmed layer, a pressure drop occurs as the atmosphere attempts to establish 
hydrostatic equilibrium, Once the surface pressure gradient is established, 


a flow from over the cool water to the warm land begins. 
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Aloft, the situation is reversed. The heating of the layer 
near the ground causes convection over the land. At higher levels in 
the atmosphere, where a stable stratification exists, the convection 
will result in a slight cooling of a deep layer aloft. The cooling 
is so weak that the increase in density is not able to cancel the 
effect of decreased density at the surface, This also,resultsiin a 
pressure gradient but opposite in sense to that near the surface and 


a flow from land to sea occurs. 


The circulation is completed with subsidence over the sea 


as a result of conservation of mass. 


Detailed observations of the sea-breeze are very scarce and, 
as a result, only the broadscale features of this circulation are 
known. Hence, it is difficult to judge the performance of a numerical 


model in its simulation of this phenomenon. 


1.3 The Calm Gradient Wind Case 

The implication of calm gradient wind is that there is no 
atmospheric motion in areas not influenced by the sea-breeze. Conrad 
(1928), on five calm summer days, noticed that the sea-breeze began over 
the coast Getween 7JOC_ALM. and 6:00 ALM. and by watching ruffling of 
the water surface, he observed its seaward progression until the sails 
OL a boat, fone, Ki lome ver Off snore, were aiteoquéed atl 9:OG;mem, Thus, 
the sea-breeze circulation spreads slowly landward and seaward during 


the day. 
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Wexler (1946) cites several examples in the tropics where 
the landward penetration of the sea-breeze reached over 150 kilometers 
from shore and was detectable 125 kilometers out to sea. He states, 
however, that in temperate latitudes the landward penetration is 
limited to 30 or 40 kilometers. Hills and forests appear to obstruct 
the movement inland. The sea-breeze of Massachusetts is retarded by 
Blue Hill, 180 meters high and 12 kilometers from the sea, while farther 


north its penetration is sreater. 


The vertical structure of the sea-—breeze circulation in the 
tropics from observations by van Bemmelen (1922) under calm gradient 


wind conditions is given in Figure (1-1). 
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FIGURE (1-1) Vertical structure of the sea~breeze at Batavia. 
Velocities are given in m/sec. 
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Since the vertical structure depends on the solar heating, cloudiness, 
latitude, and season affect the results of Figure (1-1) by changing 


the vertical extent of the circulation and the strength of the winds. 


As the day progresses, the sea-breeze and upper land-breeze 
are deflected by an ever-increasing amount from the action of the 
Coriolis force. The hourly variation of the wind during a sea-breeze 
situation at Hoek van Holland, as given by Bleeker and Schmidt (1947), 


is presented in Figure (1-2). 
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FIGURE (1-2) Hourly variation of the sea-breeze at a coastal station. 
(Hodograph) 
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1.4 The Sea-breeze with an Offshore Gradient Wind 

The sea-breeze circulation is greatly modified by the presence 
of an offshore gradient wind. A sea-breeze which develops in opposition 
to the gradient wind is termed a "Cold Front-like" sea-breeze by Defant 
(1951). The development of this sea-breeze takes place several 
kilometers out at sea and then moves landward, often not reaching the 
coast until mid-afternoon. In the morning, the offshore gradient wind 
carries warmed air from the land out to sea andjas a result, the 
temperature gradient is displaced seaward causing the sea-breeze to 
develop out at sea. As the sea-breeze passes a given point on the land, 
that location experiences a marked windshift accompanied by a drop in 
temperature. Dince these are the characteristics of a-cold frontal 


passage, this form.of the sea-breeze is well named. 


The sea-breeze front does not move so far inland as the sea- 
breeze in the calm gradient wind case. In fact, the sea-breeze front 
may not even reach the land or if it does, it can be moved "backwards" 


by the opposing wind. 


Observations by Koschmieder and Hornickel (1942) at Danzig 
indicate that the vertical extent of the frontal sea-breeze averaged 
400 meters and occasionally reached ‘700 meters while, in the calm 


gradient wind case, the depth did not exceed 200 meters. 


1.5 The Onshore Gradient Wind Case 
When the gradient wind is onshore, a strong rise in temperature 


of the air over land is prevented and the sea-breeze component of the 
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wind is then very weak. Sutcliffe (1937) found statistically at 
Felixstowe, England, that with a gradient wind from the sea on clear 
days, the added sea-breeze effect from 6:00 A.M. to about 2:00 P.M. 


was approximately 1 m/sec in the lowest 300 meters. 


1.6 The Land-breeze 

As @ result of incréased stability over the land during the 
night, the land-breeze is weaker than its sea-breeze counterpart. The 
typical seaward range of the land-breeze, as determined by smoke 
observations off the coast of Hngland, is ten to fifteen kilometers 


and eight kilometers off the Baltic coast. 


The vertical extent of the land-breeze is considerably smaller 
than that of the sea-breeze. Zinner (1919) observed an average depth 
of 180 meters for the land-breeze at Burgas on the Black Sea and an 
average depth of 690 meters for the sea-breeze. The contrast between 
the sea-breeze regime and the land-breeze regime is also evident from 


Figure (1-1). 


1.7 Previous Numerical Work 

The first attempt to integrate the non-linear equations for 
the sea-breeze was that of Pearce (1955) but his model used an 
unrealistic mechanism for the input of heat to the atmosphere. In 
spite of this assumption, the sea-breeze circulation developed and 


the centre of the circulation moved inland during the day. 
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Using a sinusoidal variation with time for the land temperature 
and using diffusion to heat the atmosphere, Fisher (1961) modelled the 
sea-breeze but did not account for changes in stability in the low 
levels of the atmosphere. His results were in error after a few hours 


because of computational instability. 


Estoque (1961) published solutions of the sea-breeze with a 
model that incorporated stability changes but assumed hydrostatic 
equilibrium and only approximate mass conservation. In spite of these 
restrictions, this model gave excellent results in that it essentially 
duplicated the observations for the sea-breeze perfectly but was 
unsuecessiul in describing the land-breeze phase of the circulation. 
Estogue's model seems plagued by computational instability. The author's 
version of Estoque's model was stable only if strong smoothing was used 


and Moroz (1968) found Estoque's model to be unstable after ten hours. 


McPherson (1970) used Estoque's model in three dimensions to 
study the affect of a coastal bay on the sea-breeze. The results of 
this model showed the distribution of vertical velocity around the 
bay. McPherson admits the use of a filter to keep the model stable 


but did not publish solutions after 12 hours of meteorological time. 


The most significant contribution by Estoque is the method 
of computation near the ground and details of this are supplied in 


Chapter 2. 


The model used by Neumann and Mahrer (1971) incorporates an 


improved sinusoidal heating function for the land, does not assume 
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hydrostatic equilibrium and gives essentially perfect mass conservation. 
This model is very successful in duplicating the existing sea-breeze 
observations. In addition, it is the first model to handle the land-~ 
breeze part of the circulation adequately. It is this model that forms 


the basis for the results presented in the present study. 
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CHAPTER 2 


THE MODEL 


Ze. Introduction 

In this chapter, the equations which define the sea-breeze 
model are presented. Ultimately, this set of equations will be solved 
numerically by use of finite-differences on a network of grid points. 


Details of the grid will be given in Chapter 3. 


In devising a model to study the sea-breeze circulation, one 
is faced with the problem of the strong vertical gradients of wind and 
temperature near the earth's surface. As a result of these rapid 
vertical changes, it is necessary to increase the number of vertical 
grid points in order to obtain satisfactory resolution near the ground 


in a study of this region using a grid point method. 


Fisher (1961) placed the lower boundary of his grid at the 
earth's surface and used a large number of grid points to deal with 
the large vertical gradients of wind and temperature near the ground. 
However, Fisher made no attempt to include the effects of static 


stability. 


In contrast, Estoque (1961) placed the lower boundary of his 
grid above the region of strong gradients and used an interpolation 
scheme to compute the profiles of temperature and wind in the boundary 
layer below the grid. With this method, Estoque was able to account 


for changes in the wind and temperature profiles as well as changes in 
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friction, resulting from changes in static stability in and below the 
grid. This interpolation method will be discussed fully in later 


sections. 


In this study, Estoque's approach was used since it has the 
ability to distinguish between regimes of forced and free convection 


and because it requires fewer numerical computations. 


2.2 Partitioning the Atmosphere 

Following Estoque (1961), the region of the atmosphere that 
is influenced by the sea-breeze is divided into two sub-layers. The 
lower sub-layer extends from the earth's surface to the top of the 
constant flux layer which is assumed to be at a height of 50 meters. 
The upper sub-layer extends from 50 meters to a height where the effect 
of the sea-breeze vanishes. Estoque (1961) states that observations 
suggest a value of two kilometers for this height. For numerical 


convenience, a height of 1950 meters is used. 


Estoque (1961) further states that the horizontal extent of 
the region influenced by the sea-breeze is about one hundred times 
its height. Again, for convenience, a value of 140 kilometers is 


used. 


The upper sub-layer is divided into a network of grid points. 
The arrangement of the grid points is in equally spaced rows and 
columns with the lowest row at a height of 50 meters and the top row 
at 1950 meters. The columns are spaced so that the extreme left and 


right columns of the grid are coincident with the 1a’ al boundaries 
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of the sea-breeze region outlined in Figure (2-1). 


No influence of the sea-breeze at this height 


1950 


Upper sub-layer 


Height (m.) 


Lower sub-layer (Constant flux layer) 


FIGURE (2-1) Partitions of the atmosphere. 


2-3 The Co-ordinate System 

The co-ordinate system used is a local tangent plane 
Cartesian system. It is assumed that the coastline is infinitely long 
and runs in a north-south direction with the land to the east and the 
sea to the west. The positive "x" axis points eastward (from sea to 
land), the positive "y" axis points northward and the positive "z" 


axis points vertically upward. 


The above orientation of the coastline is chosen for 
convenience only and the sea-breezes described by the model are valid 


with any orientation of the coastline. 
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The equations that are solved numerically to predict the 
sea-breeze circulation in the upper sub-layer are basically the Navier- 
Stokes equations for an incompressible turbulent atmosphere on the 
rotating earth plus the equation of continuity and the thermodynamic 
equation. To close this system of equations, it is necessary to use 
the equation of state for an ideal gas and Poisson's equation defining 
potential temperature. Since the horizontal extent of the sea-breeze 
region is small compared.to the radius of the earth, the equations can 
be written in local tangent plane co-ordinates. Neglectimg radiation 


and atmospheric moisture, these equations can be written as follows: 


O = WGA (3.4.5) 


where V = three dimensional mean wind vector (=Lu + yar +kuy ) 


= specific volume for dry, air 


= pressure 
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eddy diffusivity for momentum 

= potential temperature 
=ospecificagas constahtrfor dny*tair 
= temperature 

angular velocity of the earth 


= acceleration due to gravity 
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reference pressure (1000 mb) 


~ 
Saas 
et 


eddy diffusivity for heat 


Zee nes ddy DitiusivitLes 
From the mixing length theory, Hess (1959) states that Prandtl 


defines the turbulent shear stress: 


and, as a result, the momentum transport by turbulence can be written: 


0 
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Sutton (1953) generalizes the above by stating that the transport of 
a ice eile property by eddies is given by: 
Gai sane 69 lesbo pe MObetfty ithe | 
thie’ 3a \K Se) & (x ls aie Si) (3.5.1) 
where \X = scalar fluid property 

Kx = eddy diffusivity in "x" direction 

Ky 

Ke 
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eddy diffusivity in "y" direction 
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eddy diffusivity in "z" direction 
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In this study, the quantities being transported by eddies are 
the horizontal momentum of the wind and the potential temperature. 
Outside the boundary layer, the changes in the horizontal gradients of 
these quantities are far less than those in the vertical and as a result, 


only the vertical term of (2.5.1) appears in the prediction equations. 


From the expressions for the eddy diffusivities given by 
Lumley and Panofsky (1964), it can be seen that they are a maximum at 
the top of the constant flux layer. The magnitudes of the boundary 
layer eddy diffusivities drop to near zero at: sufficiently great heights 
in the atmosphere where ground induced turbulence vanishes. The manner 
in which the diffusivities decrease with height is not known. lstogue 
(1961) used a linear decrease and McPherson (1970) assumed an exponential 


drop. in this suucy, a linear decrease is used yielding; 


K(2) = K(s) = (.5.2) 


where K eddy diffusivity 


= 
I 


height of constant flux layer 


am 
il 


height at which eddies vanish | 


It is required that the height at which the eddies vanish is 
the same as the height at which the influence of the sea-breeze is 
negligible. Hence, H is set equal to 1950 meters and as previously 


mentioned, the top of the constant flux layer is 50 meters. 


Iumley and Panofsky (1964) state that over a large range of 
stabilities the diffusivity for heat, Ky » and the diffusivity for 
momentum, Ken , are not significantly different from each other. However, 


they also point out that this question is highly controversial. For the 
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lack of a better assumption, therefore, it will be assumed that Ky 


and Ky, are equal, 


2.6 The Assumption of Incompressibility 
Batchelor (1967) derives the following three conditions 


necessary for air to behave as if incompressible: 


% 


= hee) Gury 
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where = (fatlo of theyspecitic heats (for air 


the dominant frequency of oscillations in the flow field 


iH 


speed of sound 


HT 


pressure 


density 
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The spatial distributions of the wind vector U are characterized 
by a length scale L (meaning that in general U varies only slightly over 
distances small compared with L) and that the variations of the magnitude 


of U with respect to position and time have the magnitude U. 


Equation (2.6.1) states that the flow velocity must be small 


compared to the speed of sound. The second equation (2.6.2) states the 
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frequency of any oscillatory motion in the air must be small compared 
to the frequency of sound and (2.6.3) restricts the vertical scale of 
motion to be much less than the "scale height" (defined as ay ) 
of the atmosphere. It will be shown a posteriori that the model, 


outlined in this chapter, satisfies the above conditions. 


2e/ The Lower Sub-layer 


In the lower sub-layer, the vertical fluxes of heat and 
momentum are assumed to be constant in the vertical. Hence: 
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potential temperature 


The purpose of the lower sub-layer equations is to provide 
the lower boundary conditions (necessary for computation of the 
vertical derivatives of the potential temperature and wind) for use 


by the prediction equations of the upper sub-layer. 


In practicey the equations for the lower sub-hayer/ are used 
to "adjust" the values of temperature and wind at the lower boundary 
of the grid (z=h) to ensure that these values are consistent with the 
profiles, the values present at the ground (z=0) and the values present 


at the second row of grid points (z=h+Az where Az is the spacing 
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between the grid rows). 


In order for the prediction equations in the upper sub-layer 
to use the values produced by the boundary layer equations, continuity 
in the wind and temperature and their vertical first derivatives is 
required. In order to accomplish this, the constancy of the heat and 
momentum flux is extended from the top of the lower sub-layer to the 
first interior grid row of the upper sub-layer, that is, to a height 


of ht+tAz. 


Since both statically stable and unstable regimes may exist 
in the lower sub-layer, two sets of boundary layer equations are 
required. The gradient Richardson number is used as the indicator of 
stability. The eatin fron stable to unstable and from unstable 
to stable is assumed to occur at a critical Richardson number Ric : 


The unstable regime is assumed to occur with: 


Ri = Ri. Coe) 


and the stable regime with: 


PEL ene (3.4) 


where the Richardson number is evaluated at a=(h+A z)/2 and defined as: 
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Following Priestley (1959), the critical value of the 


Richardson number is taken to be -0.03. 


2.8 The Unstable Regime 

Using dimensional analysis and the assumption of a constant 
heat flux, Lumley and Panofsky derive the following equations for 
determining the temperature profile and eddy diffusivity in an unstable 


boundary layer: 
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where ) = average potential temperature in the constant flux layer 
@ =a non-dimensional constant 

B = an integration constant 

Q = air density 

Gg = acceleration due to gravity 

H = heat flux 


Cp = specific heat of air at constant pressure 


K = eddy diffusivity 


Starting from these equations, Mahrer (personal communication) 
derives an expression for the eddy diffusivity at the top of the lower 


sub-layer in terms of the values at the ground and those at the second 
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Although Priestley (1959) states that (2.8.3) is not valid 
below one meter, it will be assumed to be valid at z=z, where z, is 
the roughness height. The potential temperature at the roughness 
height is set equal to the surface temperature. In addition, since 
Lumley and Panofsky state that the value of the constant C is near 


unity, it was set equal to one. 


From (2.8.3) it can be written that: 
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where the subscript h+Az refers to values at the second grid row and 
the subscript o refers to values at the ground. Subtraction of (Orava) 


and (2.8.5) gives: 
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and substitution into (2.8.3) yields the diffusivity at z=h: 
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To determine the values of U,V, and 6 at the lower boundary 


of the upper sub-layer (z=h), the one-third power law of (2.8.2) is used 


in the. £orm : 
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following result for the potential temperature at z=h: 
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(3.3.9) 
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If it is assumed the profiles of wind and temperature are similar, then: 
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since the wind at z=z, is zero. 
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The components of the wind follow 
logically from (2.8.10): 
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A} The Stable Regime 
Estoque (1959) derived the following expressions for the wind, 


temperature and diffusivity at the lower boundary of the upper sub-layer: 


a. 
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i ame las 5 it (3.9.2) 
OMe cent (REA Ma (3.4.3) 
where K, = von Karman's constant = 0.4 


Qu = Richardson number 


A = 2 constant = —0303 
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The values of the W and ‘VU components at z=h are computed from (2.8.11) 


end (2.6312 )s 


There is a marked discontinuity in the results for the stable 
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2.10 The Boundary Conditions 
Following Neumann and Mahrer (1971), the following boundary 


conditions are used. Since the sea-breeze has no influence on the top 
boundary of the upper sub-layer, the meteorological variables will not 


change with time at this point. As a result for the top boundary (z=H): 


2 (u,v, p, 8) 
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On the lateral boundaries, the conditions that exist at a 
large distance from the coastline (where the influence of the sea-breeze 


is negligible) are approximated. At the lateral boundaries: 
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On the lower boundary, the conditions are generated by the 
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equations for the lower sub-layer discussed previously, using the 


boundary conditions at the earth's surface. At z=Z,: 


Pes Wy = 0 (3.10.5) 
Te = sea temperature 
T(t) = land temperature 


me (TQ)T,) = coast temperature 


The coast is given the average of the land and sea temperatures to 


better represent Hsu's (1967) observations. 
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harmonics of Kuo's (1968) results for temperature change at five 


millimeters below the land surface on a clear summer's day. Kuo's 
results were scaled to give a peak-to-peak temperature change of about 


22 C. The surface temperature of the land as a function of time is: 
T(t) = To +1d.aim (15+ - 110") +3.5 aim ( 30¢ +75°) 


+O.5 aim (454 +669) +O.6 aim (6Ot- 115) (ator) 


where t is the time measured in hours from the time of equal land and 


sea temperatures (= 8:00 A.M.). 


time (hours) 


FIGURE (2-2) Form of the temperature wave applied to the land surface. 
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CHAPTER 3 


THE COMPUTER MODEL 


3.1 Intiweduction 
In the previous chapter, the theoretical model for the sea- 
breeze circulation was presented. In this chapter, the method used 


for solution by computer of the equations of Chapter 2 is outlined. 


Qae ine Predietion Touations 

IMmomder toneeduce the numgermot sieniticant figures 
required during numerical computation, one can express the horizontal 
velocity components. and potential temperature as the sum of an initial 


value and a "small" perturbation, i.e.: 
Ciee Oh. FU (3.2.1) 
/ 
UO SE PAUZIE a (Area) 


G, +O (37933) 
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where the primed quantities are the perturbations and the subscripted 


quantities refer to the initial state. 


Further, it is convenient to split the pressure into three 


parts as follows: 
) | po 
pee spe ip ep (3.2.4) 
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hydrostatic perturbation pressure 


= non-hydrostatic perturbation pressure 


0% 0 


The sum (p +p’) satisfies the hydrostatic equation: 


& (pe +p’) = 5g (ya. 5) 


The initial horizontal wind field is assumed to be geostrophic 
and as a result, the initial horizontal pressure gradient terms can be 


approximated by: 
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The initial temperature field. is a function of height only 


and is of the form: 
ee Ee aie ( 3.2.3) 


where loz = temperature at height z 


ex = initial temperature at z=0 


initial lapse rate (the initial rate of decrease of the 
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If it is further assumed that the initial horizontal wind 
field is a) funetionor hétent only, the prediction equations of 
Section 2.4 become, upon incorporation of the assumptions stated above: 
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= specific gas constant for dry air 


Coriolis parameter 
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total pressure 


The total pressure above is defined as: 


ile +p+ p (3.3.13) 


Since it has been assumed that the initial pressure may be a function 
of all three space variables, it is convenient to define an average 
initial pressure which is a function of height only, using the 


hydrostatic equation for a constant lapse rate atmosphere; 


Bo(2) = palo) += 


where Polz) = average initial pressure at height z 


(3.3.14) 


Polo)= initial pressure at height z=, 


The average pressure is used only to compute the specific volume in the 


prediction equations and asa result, the error introduced is less than 


one percent. 


Bao ne Grid 

As mentioned previously, the upper sub-layer is covered with 
a network of grid points in the form of equally spaced rows and columns. 
The lowest grid row is coincident with the lower boundary of the upper 
sub-layer (at z=h). The highest row of the grid is located at 2-H 
where the influence of the sea-breeze is assumed to be negligible. The 


first and last columns of the grid are placed at the lateral boundaries 
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of the sea-breeze region described previously. 


The grid rows are numbered consecutively upward, starting 
from one at the lowest row to M at the top row, where M is the total 
number of rows. Similarly, the columns are numbered consecutively, 
starting at one at seaward lateral boundary to L at the landward lateral 


boundary, where L is the total number of columns. 


The notation used to refer to the value of a variable at 
ACS Ven 21d “POrNt 1s x sage where the grid point in question is 


located at the intersection of the fitch column and the 4, - LOW. 


To refer to the surrounding points of a given grid point, 


a) , the notation as outlined in Figure (3-1) is used. 


FIGURE (3-1) Indexing of the grid. 


The equal spacing between the rows is Az and that between 
the columns is An where An and Az refer to distances measured in 


the real atmosphere. 
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3.4. Computation of the Hydrostatic Pressure Perturbation Field 


If the initial pressure field is assumed to be hydrostatic, 
then from (3.2.5) the following can be written : 


(Oly ih) gregiay pgs on Byeehy’ 
a = a (Bf P| (3.4.1) 


The derivative is approximated by a one-sided finite-difference scheme: 
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This equation is used in succession on each column of the 


grid where the equation is integrated downward, using 54 equal 


tO Zero on the. top row. 
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3.5 The Finite-difference Equations for vw and 0 


The finite-difference equations for “U and 0 are 


written with forward differences for the time derivatives: 
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The advection terms use "upstream" differencing: 
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~ (Us +%) ou = 7 ae (Us+U; 4) GA - Wi5) (3.5.3) 


for (uw) greater than or equal to zero and: 


“(ane’) BE * Ag (beta) (Ving ha) G54 


for (U,+W) less than zero. Similar expressions are used to replace 
7 


the advection term in the 6 equation. 

According to Haltiner (1971), the upstream difference 
approximation is computationally stable for a linear equation containing 
only whis term af : 


OE. AX (3.5.5) 
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determine the stability criteria for a complete set of non-linear 
equations; hence, it will be assumed that the above linear stability 


criteria provide a close stability approximation when applied to the 


prediction equations. 


The convection terms use "upstream" differencing as well 
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From Haltiner, the linear stability criterion for the above is: 
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By use of the preceding approximations, the finite-difference 
prediction equations for the ‘VU -component of the wind and the 


potential temperature become: 
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where 99; = initial potential temperature of the j th row. 


The above equations assume (ug+ Wis) and tn are both 
negative. If either or both of the above are positive or zero, it is 
necessary to change the corresponding advection term according to 


(3.5.3) and/or (3.5.6). 


3.6 Chorin's Method for the U and Wr Equations 
The method originated by Chorin (1968) is used to predict 


the tu and uy fields. His method allows the equations of motion 
to be used in component form instead of resorting to the vorticity- 


stream function method. 


In Chorin's scheme, an auxiliary velocity field is defined 
by removing the non-hydrostatic pressure terms from the uw’ and 
Uy equations and writing them in finite-difference form similar to 
the WJ and ; equations previously described. Approximating the 
pressure gradient term in the u! equation by a centred difference, 


the auxiliary velocities are; 
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The first pair of the above equations is divided by ee and the 
second pair by OAZ . After subtracting and making the approximation 


that : 
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where Aas = 


The left side of the above equation is the divergence of the wind field 
at the (n+1)th time step which, in accordance with the equation of 
continuity, is assumed to be zero. The first two terms of the right- 
hand side form the divergence of the auxiliary velocity field. The 
remaining two terms are a finite-difference Laplacian whose nodes are 


separated by two grid lengths. By defining: 
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equation (3.6.7) can be written: 
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Chorin solves (3.6.9) for the pressure field, using an 
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iteration procedure and by setting W on the lateral boundaries 


equal to wn" and setting Us ow On ene top and bottom boundaries 


ane 


equal to UWr™™"' in (3.6.8). Hence, it is required that W and 


nei 


us are known on the boundaries specified above. In this work, 
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a sequential relaxation method is used, setting P =0 at the 


boundaries for all time steps. 


3./ The Relaxation Procedure 

Equation (3.6.9) is solved by iteration. At each step in 
the iteration procedure, (3.6.9) is evaluated with the pressures from 
the previous step and in the case of the first iteration, step,. the 
values calculated at the previous time step are used. In general, 
these will not satisfy (3.6.9) and a residual pressure term will remain 


defined bys 
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MW. = iteration index 


The value of the pressure at the point Ceres) is adjusted so that 


the residual vanishes. This new pressure is the (m+1)th iterate: 
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Subtraction of (3.7.2) from (3.7.1) gives the iteration equation: 
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This equation can only be used in the area of the grid which is two or 


more columns or rows from the boundaries. To derive the iteration 
equation for the region one grid length away from the boundary, we 
make use of the fact that the appropriate component of the auxiliary 
velocity on the boundaries stated in Section (3.7) is replaced by the 
appropriate component of the "n+l" velocity field, which forces the 
normal derivative of the pressure across the boundaries equal to zero. 
Hence, every pressure at a point one grid length from a boundary has 


an equal image pressure lying outside the boundary. 
On Figure (3-2), the following pressures are equal: 
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FIGURE (3-2) The image pressures. 


Using the above, (3.7.1) and (3.7.3) can be used one grid length from 


the boundary. The iteration equation for the column 4 =2 is: 
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Similar expressions are used at other grid points lying one grid length 


from the boundary. 


There are four points in the corners of the above region which 


require a third iteration equation derived also by the use of images; 
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A similar expression is used for the other three corners. Figure (3-3) 
shows a corner of the grid and indicates which iteration equation is 


Used for cach of the grid points. 


Boundary 


Boundary 
FIGURE (3-3) Use of the three iteration equations. 


The iteration procedure is applied sequentially to each point of the 
grid as many times as necessary to reduce the residual at every point 
to less than some pre-determined maximum error. In practice, the 


iteration procedure converges more rapidly if an over-relaxation is 


BO 


used at each grid point. If ‘such a procedure is used, (3.7.3), (3.7.5) 


and (3.7.7) respectively become: 
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where IN = over-relaxation factor 


lie wvalue7or IX which gives the most rapid convergence is most accurately 
determined by experiment and was found to be 1.7. 
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Once the values of P have been determined, (3.6.3) to 
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27 rl n 
(3.6.6) are used to find U and US , 


3.8 Computation of the Vertical Velocity at the Lower Sub-layer Boundary 
Figure (3-4) shows a portion of the grid near the lower 


boundary. 


The equation of continuity written in centered differences 


for the point P.is: 
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FIGURE (3-4) Boundary value of the vertical velocity. 


Since the vertical velocity is zero at the earth's surface. If a linear 


increase in the wind is assumed in lower sub-layer, then: 
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Hyen if Chis ascimption ts °in error, the vertical velocity field is 
hardly affected since it is not too sensitive to the values of Wj \ 
Substitution of the above into (3.8.1) gives the result: 
n+\ h 7 2 n+ 
macaw I 7 AAs. ( De get =e) (3.9.4) 
3.9 Filtering the Velocity Fields 


A one dimensional three-point filter is applied at each time 


step to all three components of the velocity field. The form of the 
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where Fy = unfiltered values 
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a constant which determines the response of the filter 


To control non-linear instability, it is necessary to remove 
two to three gridlength waves with the filter. According to Haltiner 
(1971), the response of (3.9.1) with S=0.5 is: 

R = sayin eves: 
: Ms: 
where R =| GOspolce 
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wavelength of a wave in units of Ax 
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The response is the amplitude of a given wave after application of the 
filter divided by its amplitude before filtering. Clearly for a two 


gridlength wave, the response is zero. Figure (3-5) shows the response 


1.0 
» 75 
w 
6 
© 50 
(79) 
® 
© 25 
fe ites TE AT Te Lee Ree es es 


és fs) “} 5 6 7 8 Sheet) 
Wavelength in gridlengths 
FIGURE (3-5) Filter response. 
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for other wavelengths. 


3.10 Ordering of the Steps 
To help clarify what has been presented so far, it is helpful 
to show the steps followed in the computer model. 
a) determine initial values 
b) compute Q ; Ir aint Ww at z=h using the lower 
sub-layer equations of Sections (2.8) and (2.9) 
c= Priter ue ; ay and Uy fields using (3.9.1) 


ad) apply 2 20 boundary condition to WU’ 3 Ay and 


8 fields 

e) compute temperature using (2.4.5) 
f) compute hydrostatic pressure perturbation from (3.4.3) 
2) compute U"’™ and UW" Fields using (3.6.1) and (3.6.2) 
h) compute ste and WS from (3.5.11) and (375.12) 
i) apply temperature wave to land surface (2.10.7) 
3) SRS WN PC HeR. 
k) determine non-hydrostatic pressure perturbation using 

equations of Section (3.7) 
1) compute eee es a using (3.6.3) and (3.6.5) 
m) repeat steps b) to 1) until desired length of forecast 


is achieved 


n) display results 


A fully documented listing of the computer program is included in 


Appendix A. 
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CHAPTER 4 


THE RESULTS OF AN INTHGATION WITH THE 


ATMOSPHERE INITIALLY AT REST - THE COMPARISON MODEL 


ele ncrouuction 

The computer program, used by Neumann and Mahrer (1971), 
was modified by including the vertical Coriolis term, changing some of 
the lower sub-layer equations and increasing the peak-to-peak amplitude 


of the temperature wave applied to the land from 20 GC to 25 C. 


The purpose of duplicating Neumann's and Mahrer's results 
is to provide a standard by which to compare the results that are 
presented in succeeding chapters and for this reason, the model will 


be referred to as the "Comparison Model", 


Lee initial Conditions and Fields 

The model was integrated for a period of twelve hours 
Scartaney at o:OONkM., using a: 20urow by 29 colummegriad.. Im the 
computer program, an extra row is added to the grid to store the data 
for z=0, making the grid 21 rows by 29 columns. The finite-difference 
data concerning the grid as well as the initial lapse rate, horizontal 
wind components, surface temperature and surface pressure are presented 
in Table (4-1). From these data, the initial fields of Table (4-2) 


were computed. 


The size of the time step was chosen to satisfy (3.5.5), 
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(3.5.8) and (3.5.10). The largest expected values of the 


Uy component and eddy diffusivity were: 
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Hence, the restriction on the time step for each of these values was: 
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To have a margin of safety, the time step used was 150 seconds. 


NUMBER OF HORIZONTAL GRID POINTS 29 
NUMBER OF VERTICAL GRID POINTS 21 
HORIZONTAL SPACE INCREMENT (m) 5000.0 
VERTICAL SPACE INCREMENT (m) 100.0 
TIME STEP (sec) 150.0 
LAPSE RATE (deg/m) 0.0065 
U COMPONENT (m/sec) 0.0 
V COMPONENT (m/sec) 0.0 
SURFACE TEMPERATURE (K) 299.0 
SURFACE PRESSURE (mb) 1000.0 

TABLE (4-1) Grid data and initial conditions. 
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TABLE (4-2) 


TEMPERATURE 


(K) 


299.0 
298.7 
298.0 
297 4 
296.7 
296.1 
295.4, 
294.48 
29he1 
29365 
292.8 
22202 
291.5 
290.9 
mo0ee 
289.6 
288.9 
288.3 
287.6 
287.0 
286. 3 


POTENTIAL 
TEMPERATURE (K) 
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PRESSURE 


(mb ) 


1000.0 
9943 
983.0 
971.68 
960.7 
ood 
938.8 
928.0 
91763 
906.7 
896.2 
885.8 
875.5 
865.3 
$55.42 
845 «2 
835.42 
825 4 
815.7 
806.0 
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HEIGHT 
(m) 

0.0 
50.0 
150.0 
250.0 
350.0 
450.0 
gio ih a8) 
650.0 
750.0 
850.0 
950.0 
1050.0 
1150.0 
1250.0 
1550-20 
1450.0 
150.0 
1650.0 
1750.0 
1850.0 


1950.0 


45 


7” 


es 


cae 


- 
7 


2 oo 
7 


* 


. 


is 


7 
: _ 
“ . 
r 
7 
i - - ; 
x _— i 7 


a 


& 
OO 


. 
_ 
a 

_ 


a 
arom eee ) thd 


= 
= 
— 


[HR eR 2 
ey 


an 


46 


A latitude of 43 N was chosen for the integration since this 
-+ ~1 
allowed the Coriolis parameter to be set equal to |O sec ., Because 
the sea-breeze does not depend strongly on the latitude, this choice 


of latitude does not severely restrict the applicability of the results. 


4. 3@ Results 
The computer program displayed the wind and temperature fields 


every two hours of meteorological time, that is, after 48 time steps. 


The results of the model at 10:00 A.M. showed that the centre 
of a weak negative circulation (negative in the mathematical sense) was 
evident at a height of about 350 meters above the coastline. At this 
time, the strongest onshore wind of about 0.5 m/sec was at the coast 
and a weaker flow of larger vertical extent was present aloft. The 
turning effect of the Coriolis force was beginning to be felt since a 


small Ww -component was present. 


The noon results (Figure (4-1)), as expected, showed that 
the circulation had strengthened. The centre of the circulation, 
located 3 km inland was at a height of nearly 500 meters. As was 
previously the case, a weaker return flow existed aloft. The maximum 
ascent velocity was 3 cm/sec with the maximum descent about the same. 
The action of the Coriolis force continued to strengthen the VW -com- 
ponent of the wind and as well, a striking vertical temperature 


gradient was present near the ground. 


At 2:00 P.M., which was shortly after the time of maximum 


temperature of land surface, the results (Figure (4-2)) showed that 
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FIGURE (4-1) Results for noon, Horizontal velocities are given 
in m/sec, vertical velocities in cm/sec and temperatures in C. 
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FICURE (4~2) Results for 2:00 P.M. Horizontal velocities are given 
in m/sec, vertical velocities in en/sec and temperatures in C. 


AD 


the onshore flow had reached nearly 3 m/sec. The centre of the 
circulation was now 8 km inland at a height of about 600 m. At this 
time, the maximum descent was slightly over 4 cm/sec with the maximum 
ascent about 50 percent larger. Because of the distorted scale of the 
diagrams, continuity may not appear to be satisfied, i.e. it is not 
necessary to have large vertical velocities in order to have large 
horizontal velocities. The temperature field showed a slight upward 


bulging of the isotherms near the region of greatest upward motion. 


Although the land had been cooling for nearly three hours, 
the results for 4:00 P.M. (Figure (4-3)) showed that the sea-breeze 
continued to intensify. The strongest onshore wind was about 4.5 m/sec. 
By this time, the Coriclieforcevhad turned’ theronshore flow nearly 15 
desrees=i ron a line perpendicular to the coast A The veitect of the 
high turbulent friction over the land was clearly shown by the strong 
horizontal lgradients of the W =component of the wind. This effect 
via continuity was responsible’ for increased ascent over the land. 

The centre of the circulation which continued to accelerate inland was 
located 20 kilometers from the coast at a height of 800 meters. The 
bulging of the isotherms in the region of strongest ascent continued 
to increase although the vertical temperature gradient near the ground 


had weakened with the continued cooling of the land. 


Figure (4-4) gives the results for 6:00 P.M. which showed 
that the sea-breeze had started to decrease in strength. The wind had 
been turned by 30 degrees from the action of the Coriolis force. The 


enhancement of the ascent by friction was even more striking than in 
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FIGURE (4-3) Results for 4:00 P.M. Horizontal velocities are given 
inim/ sec, vertical velocities in em/sec and temperatures in C. 
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FIGURE (4-4) Results for 6:00 P.M. Horizontal velocities are given 
im m/sec, vertical velocities in cm/sec and temperatures in C. 
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FIGURE (4-5) Results for 8:00 P.M, 
in m/sec, vertical velocities in cm/sec and temperatures in C. 
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the previous results. There had been so much cooling over the land 
that an inversion was present. Although the horizontal and vertical 
velocities had decreased only slightly, the core of the circulation 


was 30 kilometers inland at a height of approximately 900 meters. 


The results for 8:00 P.M. (Figure (4-5)) showed that the sea- 
breeze was dissipating with an onshore flow of nearly4.5 m/sec, 
turned 40 degrees by the action of the Coriolis force. The vertical 
velocities are probably in error since) the core was near the lateral 
boundary. The temperature field showed a strengthened inversion over 


land and an isothermal layer above it to a height of 400 to 500 meters. 


Figure (4-6) gives the hourly variation of the wind at the 


coastline for the twelve hour period of integration (cf. Figure (1-2)). 
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FIGURE (4-6) Predicted hourly variation of the wind. 
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4.4 Verification of the Incompressibility Assumption 
Batchelor's criteria were applied to the results given by 


the Comparison Model. Condition (2.6.1) states that: 
> 
<a 


ic 


Since the strongest wind was about 5 m/sec, (2.6.1) is verified. 


Since the smallest time scale that could be resolved in this model is: 


VAL = 300 sec * = 


n 


the inequality (2.6.2) is satisfied. The most stringent condition is 
(256.3): 


egh area (4.4.1) 


Biers tiGetenecesssry. vo derermine the scale length Lb’: 


ee 
lee ae 
Oz 
Evaluating (4.4.1) at 4:00 P.M. gives: 
ee TUCO eR 


Hence, (2.6.3) becomes: 


ggk _ @nO.ba 4 
Bp 


From the application of Batchelor's criteria, it can be seen that the 


assumption of incompressibility is reasonable. 
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A check was made to see whether the model conserves energy. 


The first law of thermodynamics in the form: 


ais 6, lis See ‘2 
re an or (4.5.1) 


where Q = heat input 
‘| = temperature 
Cp = specific heat at constant pressure 


A 


specific volume 


Pp = pressure 


| 


= time 


was integrated over the entire grid with the exception of the boundary 


points. Equation (4.5.1) can be rewritten in terms of potential 


temperature as: 


ah We eile, 
a = chal it Gas 20) 


and upon substitution of (2.4.3) becomes: 


ene t eo cle 


Substitution of (4.5.3) into (4.5.1) and expansion,of the total 
derivatives gives; 
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ee (4.5.4) integrated over the grid yields: 


ths tom D 
5 CP a ( K wa Dec AIS = 


i= oe 
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“I (5) 
x wae pes \ ee ar Neen 
AZO A=3 L=2 ae 
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=a A=3 AeA AO 
(4.5.5) 


where /\t = time step 


AX = horizontal space increment 


Az 
M 


vertical space increment 


i 


number of grid rows 


L. = number of grid columns 


Term (1) of the above is the heat entering the grid through 
- the lower boundary by turbulent diffusion, term (2) is the increase 
in enthalpy, term (3) is the heat leaving the grid through the 
lateral boundaries, term (4) is the heat leaving the lower boundary 
of the grid by convection, term (5) is the change in potential energy 
in the pressure field and term (6) is the production of mean and 


turbulent kinetic energy via the action of the pressure field. The 


(22+). 
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major part of term (7) is the increase in potential energy and the 


rest of this term is the production of kinetic energy. 


Figure (4-7) gives the cumlative magnitude of each of the 
terms of (4.5.5) since $:00 A.M. Terms (5) and (6) are not included 
since they are small. The maximum imbalance of equation (4.5.5) 
occurred at the twelfth hour when the right hand side was six percent 


larger than the left hand side. 


i 
The centre of the circulation is located at the intersection 
of the zero isopleth of the u-component and the zero isopleth of the 
vertical velocity. 
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FIGURE (4-7) Cumulative magnitydes oh the terms of the energy equation. 


The units of energy are 10 Jm kg 


CHAPTER 5 


SYNOPTIC CONTROL OF THE SEA-BREEZE 


5.1 Introduction 

The motionless atmosphere in which the sea-breeze of the 
previous chapter occurred is rarely observed. In an attempt to 
approach physical reality more closely, the effects of initial 
atmospheric motion, initial lapse rate and ground heating on the sea- 


breeze circulation were studied. 


Estoque (1962) presented the results of similar solutions 
using a model which assumed hydrostatic equilibrium and which also 


used only an approximation to the incompressible continuity equation. 


Bee ine Biiectsvot incressed sotabulicy 


The initial conditions and grid parameters used for this 
integration are given in Table (5-1). These data were the same as 
those for the Comparison Model except that the initial lapse rate had 
been made more stable. From the initial conditions, the initial fields 
of Table (5-2) were computed. To prevent an inordinate number of 


figures, only the results for 10:00 A.M. and 4:00 P.M. are presented. 


The results of this solution parallel those of the previous 
chapter except for a few interesting details. The 10:00 A.M. results 
(Figure (5-1)) were rather unexpected since they showed that the 


increased stability had strengthened both the onshore wind and the 
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NUMBER OF HORIZONTAL GRID POINTS 
NUMBER OF VERTICAL GRID POINTS 
HORIZONTAL SPACE INCREMENT (m) 
VERTICAL SPACE INCREMENT (m) 
TIME STEP (sec) 

LAPSE RATE (deg/m) 

U COMPONENT (m/sec) 

V COMPONENT (m/sec) 

SURFACE TEMPERATURE (K) 


SURFACE PRESSURE (mb) 


TABLE (5-1) Initial conditions and grid parameters for the stable 


lapse rate case. 


return flow aloft. The stable air above the unstable layer near the 


ground produced a sharp vertical change in the lapse rate in the low 


2°6 


levels. Due to the Cae 


relatively strong rise in temperature occurred, which resulted in an 
increased onshore flow. Aloft, the stable stratification produced 
strong cooling thus causing an increased return flow as well. 


the vertical penetration of the landward flow had been reduced. 


In contrast, the 4:00 P.M. results (Figure (5-2)) showed 
that the intensity of the sea-breeze had been diminished by the 
increased initial stability. The centre of the circulation was 
located 9 km from the coast at a height of approximately 500 meters. 


The strength of the sea-breeze and vertical velocity were less while 
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ROW TEMP ERATURE POTENTIAL PRESSURE HEIGHT 


(K) TEMPERATURE (mb ) (m) 
(K) 
a4 299.0 299.0 1000.0 0.0 
a 298.9 EAS be Ia 994-3 50.0 
3 298.8 300.3 983.0 150.0 
Js 298.7 301.2 971.8 250.0 
5 298.6 302.1 960.8 350.0 
6 298.5 303.0 949.9 450.0 
he 298.4, 303.9 939.1 550.0 
8 298.3 304.8 928.4 650.0 
9 298.2 305.7 917.8 750.0 
10 298.1 306.6 907.4 850.0 
all 298.0 BUVeo 897.0 950.0 
1 297.9 308.4 886.8 1050.0 
ie 297.8 309.3 876.7 1150.0 
14 29 Tek 310.2 866.7 1250.0 
15 297.6 oa laa 856.8 1350.0 
16 297.5 312),0 84.7.0 1450.0 
17 297.4 312.9 837.4 1550.0 
18 297 3 313.9 827.8 1650.0 
19 297.2 314.8 818.4 1750.0 
20 27 (ou 315.7 809.0 1850.0 
pal 297.0 316.6 799.8 1950.0 


TABLE (5-2) Initial fields for the stable lapse rate case. 
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10:00 A.M. results for the stable case. Horizontal velocities 
are given in m/sec, vertical velocities in cm/sec and temperatures in C. 
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FIGURE (5-2) 4:00 P.M. results for the stable case. Horizontal velocities 
are given in m/sec, vertical velocities in cm/sec and temperatures in C. 
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the return flow aloft was greater than in the results for a less stable 
atmosphere. The flow was less than in the Comparison Model because 

the rapid vertical change in the lapse rate had been removed by 
diffusion. In addition, the heat input to the atmosphere was used to 
overcome the increased density aloft, produced by the strong cooling 
which resulted from the stable lapse rate. This effect via the hydro- 


static equation caused a relatively weak onshore flow. 


Deo AndOUEshore Plow of 3 m/sec 


NUMBER OF HORIZONTAL GRID POINTS 29 
NUMBER OF VERTICAL GRID POINTS Zu 
HORIZONTAL SPACE INCREMENT (m) 5000.0 
VERTICAL SPACE INCREMENT (m) 100.0 
TIME STEP (sec) 150.0 
LAPSE RATE (deg/m) 0.0065 
U COMPONENT (m/sec) Se 
V COMPONENT (m/sec) 0.0 
SURFACE TEMPERATURE (K) 299.0 
SURFACE PRESSURE (mb) 1000.0 


TABLE (5-3) Initial conditions and grid data for the offshore gradient 
wind case. 


Table (5-3) gives the grid parameters and initial conditions 
used for this solution. The initial fields were the same as those 


given in Table (4-2). Only the results for 10:00 A.M. and 4:00 P.M. 


are included. 
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IGURE (5-3) 10:00 A.M. results for an offshore wind. Horizontal velocities 
re. glven in m/sec, vertical velocities in cm/sec and temperatures in C. 
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GURE (5-4) 4:00 P.M. results for an offshore wind. Horizontal velocities 
e given in m/sec, vertical velocities in cm/sec and temperatures in C. 
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The results for 10:00 A.M. anaes (5-3)) showed that the 
sea-breeze had not yet been able to counteract the offshore flow. 
The vertical velocities, however, had a similar spatial distribution 
to those of the Comparison Model with the exception that the line of 


zero vertical velocity had been displaced nearly 10 km seaward. 


Shortly after 12:00 A.M., a landward flow began over the 
sea 5 km from the coast and at approximately 1:00 P.M., the sea-breeze 


reached the coast. 


Figure (5-4) gives the results for 4:00 P.M. A sea-breeze 
front, similar to that described in Chapter 1, was located nearly 10 
km inland from the coast. The combined effects of turbulent friction 
and convergence at the front produced an area of strong ascent while 
only weak subsidence occurred over the sea. As expected, the landward 
penetration, vertical extent and strength of the sea-breeze were 
retarded by the offshore flow but the ascent and return flow aloft 


were strengthened by the opposing wind. 


5 fm Ang@Onshore Flow of 3 m/sec 


With the exception of the initial UW -component, the initial 
conditions, initial fields and grid parameters were the same as the 
offshore case. The initial conditions and grid data are given in 


Tables (5—/are 


The results for 10:00 A.M. (Figure (5-5)) showed that a weak 


sea-breeze circulation had begun but the vertical velocity pattern was 


So@redaion od 


me os 
wokt cia 


=D 7 i * 
noltudingerh tabteqe waiinte e bad 


Yo satfL ocd held tipbihicne tii 2 <— dd: i is Paya esl to-6 oils ‘or 


- aa oath 


i - ; is 
cbepumon mf OL ¢iseen beosigetty Sela po “@itcoley isos: 2 
- : oe 

ei Woll truevual o ash Ousse oda wtsode 


oo 
eikery 


exeotd-noa aid Mf. OO:f yieganivoiags va fee taece oxit wont wh ae 


ae he <tuson ori} bad 


i“ 4 7 
fous A MAE O02), 792 os Liters aid! gvhy ()~2) oauye 


‘ 
OL y¥inege Hefbood iby Sf satqatS at Ladi: ends beds od saltets tag t 
} iy ; a ‘ 
A ‘> : 


notéoiet doofidtird te etedtie hetiidies art tne aiid ora? bon of 


ma z 


eltiw doesae poowds to Bek as 5 aha oped odd ts eure ed bas 
baswbasl edt .bedseqxs se /ods ter betitmse womobtedan inew Y re 


a » 
f ecosidpse ods “tai MS are: bie tnatxe Leavis leatdnciid 
: : wv 


- ‘ ‘ é ; =) 

Miele vell surest bas dpedes ‘of “ent wolt sroledic edd ‘a bebiads 
; ot 

Riango ety ot Be nedayionse as " 


+, 


oft 


wae 


¢ 
an 


: gae\e € to welt ezadpe? Sh ME 
fatsint ed) gdretiogsos- av 1 sane ~ to aot teaenco ent ce 7 
of? ae eden aut wes xatonging bry bus 2bfokt iabtint nf sets art 


Sa - va } an 
at. seviy ora web pity tne east fates tame sa. vee 
rT 


68 


1850 


1650 


1450 


1250 


1050 


850 


Height (M.) 


650 


450 


250 


50 


—40 20 0 20 40 60 


Landward Distance from the Coast (Km.) 


HORIZONTAL VELOCITY 


——U Component ---V Component 


Cee ets Wo Sea at So ee ee ee 


Height (M.) 
ao 
On 
oO 


—40 —20 0 20 40 60 


Landward Distance from the Coast (Km.) 


TEMPERATURE AND VERTICAL VELOCITY 
—— Vertical Velocity 
--— Temperature 


‘IGURE (5-5) 10:00 A.M. results for an onshore wind. Horizontal velocities 
ire -gaven ain m/sec, vertical velocities in cm/sec and temperatures in C. 
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FIGURE (5-6) 2:00 P.M. results for an onshore wind. Horizontal velocities 
are given in m/sec, vertical velocities in cm/sec and temperatures in C. 
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displaced landward by the onshore gradient wind. 


NUMBER OF HORIZONTAL GRID POINTS 29 
NUMBER OF VERTICAL GRID POINTS mal 
HORIZONTAL SPACE INCREMENT (m) 5000.0 
VERTICAL SPACE INCREMENT (m) 100.0 
TIME STEP (sec) 150.0 
LAPSE RATE (deg/m) 0.0065 
U COMPONENT (m/sec) 3.0 
V COMPONENT (m/sec) 0.0 
SURFACE TEMPERATURE (K) 299.0 
SURFACE PRESSURE (mb) 1000.0 


TABLE (5-4) Initial conditions and grid data for the onshore gradient 
wind case. 

The results for 2:00 P.M. (Figure (5-6)) showed that the 
familiar pattern of vertical velocity was displaced far inland. In 
addition, vertical velocity over the sea was nearly of the same 
magnitude as that over land. The maximum onshore flow at 2:00 P.M. 
was about 4.5 m/sec and was located over 20 km inland, After 2:00 P.M., 
the results were probably in error due to the approach of the circulation 


to the landward lateral boundary where Ur is constrained to be zero. 


5.5 The Effect of Reduced Insolation 
Since radiation is implied in this model through the 
temperature wave applied to the land, it is possible to study the 


effect of reduced heating by decreasing the amplitude of the temperature 
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TIGURE (5-7) 4:00 P.M. results for reduced heating. Horizontal velocities 
are given in m/sec, vertical velocities in cm/sec and temperatures in C. 
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The model was run with the same initial conditions and grid 
data as the integration of Chapter 4 with the exception that the 
heating of the land was halved. The results for 4:00 P.M. are given 
in Figure (5-7). Although the strength of the horizontal wind field 
was approximately half that of the full heating case, the vertical 
velocities are greatly reduced as a consequence of the much lower 


turbulent friction which has a non-linear temperature dependence. 


5.6 Summary 
Increased stability generally decreases the strength of the 


sea-breeze since stability inhibits the diffusion of heat energy into 
the atmosphere. Similarly, decreasing the ground heating reduces the 
ime ivence sot wthe sea=preezescince Jess heat cnerey is available for 
conversion into mechanical energy. The results show also that a 
linear relation between heating and stability and the strength of the 


circulation does not exist. 


The complex non-linear relation between the prevailing 
synoptic wind pattern and the sea-breeze circulation is shown by the 


results for offshore and onshore initial winds. 
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CHAPTER 6 


THE INCLUSION OF MOISTURE IN THE MODEL 


6.1 Introduction 

Since water can exist in two states at above freezing 
temperatures considered in this study, it is necessary to use two 
additional equations in the model; one for water vapour and the other 
for liquidywater.@ dn additionjeivisenecessaryitolmedify the 
thermodynamic equation to account for the transfer of latent heat 


produced during phase changes. 


6.2 The Prediction Equations for Water Vapour and Liquid Water 
The quantity used to describe the concentration of both the 


vapour and liquid is the maxing ratio which is defined as the ratio 
of the mass of liquid or vapour to the mass of dry air containing the 
vapour or liquid: 


r mass of liquid water 
ue “ mass of dry air 


. _ mass of water vapour 
a, einsss sol sdryeair 


where gy =macing-ratio-for-laqudd water 


Yay = mixing ratio for water vapour 


The rate of change of the mixing ratio with time following 
the air motion (neglecting the fall-out of water droplets since this 


is small) is caused by diffusion and condensation/evaporation, resulting 
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from vertical motion and a correction term for ensuring that the air 
does not become supersaturated or subsaturated with liquid water 
present. The correction terms result from the fact that the time 
step used is too large to adequately represent the processes of 


evaporation and condensation. Hence: 


G 
gre S Dag NAL age HC og (6.2.3) 
Dire thet WW fact Cry (6.2.4) 
Dt 
where ) = rate of change due to diffusion 


W = rate of change due to vertical velocity 


e. = correction, term 


If it is assumed that the eddy diffusivities for heat, water 


vapour and liquid kana are equal, then the diiftusion terms are: 


Dap = 2. (kK 8) KC, PSPS) 


2 (kK ete) (6 Se) 


i 


Du 


Tie meena shows that diffusion is more important than drop- 
let fall-out. Figure (6-1) shows a grid box bounded on the top by a 


mixing ratio ot 1 em/kg and on the bottom by a mixing ratio of 0 em/kg. 


= 1gm/kg 


r=0 


FIGURE (6-1) Grid box for diffusion/fall-out calculation. 


The increase of the mixing ratio at the point P from diffusion is: 
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Petterssen (1959) gives the mean fog droplet size as 10-15 p which 
have a terminal velocity of 3 mm/sec. The increase in the mixing 


ratio due to fall-out is: 
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Hence, diffusion is more important than fall-out. 


In the atmosphere, the amount of water changing phase per 
unit time following the motion is the time rate of change of the 
saturation mixing ratio for water vapour. To account for mixing ratio 
changes resulting from vertical motion, it is necessary to express the 
rate of change of the saturation mixing ratio as a function of the 


vertical velocity: 

Ae -4 (uy, T, 5)= We (6.3.7) 
where Vay, = saturation mixing ratio for water vapour 
Since the total mass of water substance must be conserved: 


Wie = 2 Whe (6.2.3) 


From Hess (1959), the saturation mixing ratio is approximately related 


to the saturation vapour pressure by: 
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ns Rs 
insg OE S ae 
where Yayg = saturation mixing ratio 
€ =ratio of molecular weight of water to the effective 


molecular weight. of dry air 


Q, = saturation vapour pressure 


— 


P 


total pressure 


Hess (1959) also gives the saturation vapour pressure as a function 


of temperature: 


y Lets al _ 
esc Gligorn =a yee (AAs) 


where Cs = saturation vapour pressure 
\_ = latent heat of condensation (assumed constant) 
Re specific gas constant for water vapour 


| = absolute temperature 


In order to begin to find an expression for (6.2.7), (6.2.9) is 


differentiated to obtain: 
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Substitution of the hydrostatic equation in the form (ignoring pressure 


changes resulting from the weight of liquid water and water vapour ) 


als a 
a x4 ig (sear ns) 
Be eh 

oie 


gives: 


exgandng wma 0 ane oltatecrt (os is te 


mate Ses * wottew ‘te, 30 haw w cata a 
“is yb 2 heated vivo oon ! — 


ott ia 


a 
Sees mewey notion - = oo) s 4 


_ . 

; awe on fateh = 4 F 
i | er apt 

t 2 em oriseery aiegevy colt niceine ads is cals 


s 


iz 
_, 
- >_* 
= 


a 
od 


. v fi | om ee, yeh - 
Cirary, | i ire! a { an Pie an 


a 


% 

& 
a 
— 


ey a ideas’ wisenmiee = 


(torcation basuane) nohieesobaed to Sand fastal= J 


J 7 ® 
i . 


Tveqev  sdew adi gaalende ep ol} foege = ie kan 
’ : i 

 suitaveumet stufoads = T a 

| 

= ' 


ek (@.8.9) (V5.0). sot Kotevens (as feeitt- a nityod” oF vabro Fl 
gett ot bedatse gi 

i a ? 

(ne ald = L. ab 2b - 
CG } Fb 22 


> 2 at 


~ 


uy 
0 a at . yents Does 
as 4 » Ml ‘ , 
} 9 a te i 


' 
7 
: _ 


ie 


Thus, (6.2.11) becomes: 


| dims y | des ye GUS (6.3.13) 
Differentiating (6.2.10) yields the Clausius-Clapeyron relation: 
Mee tcCce Se uresiam, Ac, (race) 


Gar at Ruger conte a 

and substitution into (6.2.13) produces: 
ae d tass =: male g! olf sep uted (ex f5) 
Fey ae ee lea cc RT 


The first law of thermodynamics for an ideal gas (dry air) can be 


written: 


dQ ap ee ge dt = ies ap (Crile) 


Cy = heat transferred to the gas 


Cp= specific heat capacity of dry air at constant pressure 
A = specific volume 

Assuming that the condensation takes place as a result of saturated 
adiabatic expansion, (6.2.16) becomes: 
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Flimination of between (6.2.17) and (6.2.15) gives: 
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Hence by (6.2.7) and (6.2.8) 


Vile > a= Sib (Grady 


We 2) Sys fs ANS) 
Because of the remarks prior to (6.2.17), equation (6.2.18) ic valid 
only if a phase change takes place, that is,) only if the air is 
already saturated. If the vertical motion does not result ina 


Piase change, then S=0. 


In order to: correct the errors caused by the time step 
being too large, the values of Vay and Yury at each grid point 
arewcnecked at every timevsiep to ensure that: 

a) the air is not supersaturated 


b) the air is not subsaturated with liquid water present 


If condition (a) is not satisfied, the value of {xr is lowered to 
Vass and the liquid water mixing ratio is increased by the same 
amount. Since this correction is computed over a time interval, 
“ste it is necessary to divide the above result by At im order’ to 
obtain’ the change in the mixing ratio per unit time. Hence, the 


correction term becomes for Vay >Mayg : 


Ces ieee) / At (eect) 
Chee aoe tay AE (4.9.93) 


If condition (b) is not satisfied, the air is brought to saturation 
and the mixing ratio of the liquid water is lowered by the same amount; 


however, the total correction cannot be greater than the total amount 


neta. sk (bss 8) aaa ts 
2 ths edt It Yino ot Jelld yedeip aatariapesits « a 
arcs 
i i ica a 
cote sollagh swe lane die alt Helge te econ ey 
satin rg dA a a Seeoeed in 
hast Sratsite oF qsde veld wave J 
Losmadaigmyi i et at at (a 
1a AO NE A EP (¢ 
_ oh betibhiel ot a nite sei ola 7 
eme add yd bebserott ef oltat pabite yohew Alopl edt bom ere 
Dévatict obit » pave fediqney ak Motvoerzoe atdtt ‘sole .dawoms 
ait nethao nt gh tied aahde’ vl “ad Weindingat as SEGRE 
add ovr sents Shi age Oto yale ale nt oper lt abate | 
f ohadVE@) 202 sonoced eras woltoerx0: 


2) aA\ Gat a) ae 
ah SD\ (el wt tw 8 


a ts iiitiezine con at (dy abbott ag 


79 


of liquid water present. The correction terms are for A= tars 


and Van 210: 
Cte rman \ tests 5 Me | ( b -ees) 
Gos eal A OATES AO \ east is | (&. 2. 34) 


where min = the smaller of the two arguments 


The left-hand sides of (6.2.3) and (6.2.4) are expanded in 


a way analogous to (3.2.12) to obtain the prediction equations: 
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These equations, written in finite-difference form similar to (3.5.12), 


become: 
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pare “= At aa ; 
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where YW= Way 
YO3c= Yur 


The above equations assume (Us 4yand Ue are both negative. If 


either or both of the above are positive or zero, the differencing 


scheme for the advection terms mist be changed as stated in Section (3.5). 


inpadad one unematabidi pyecciteria (3.5.5 ),.(3.569) and (3.5.10) are 


valid. 


6.3 The Boundary Conditions for the Liguid Water and Water Vapour 


For the reasons stated in Section (2.10), the boundary 


conditions are: 


a) at the top boundary 
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at the lateral boundaries 


( es GEN = 0) 


b) 
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Ox (Tf) = O 


(6.3.1) 
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c) since the total water is conserved and there is no 


advection, the boundary condition at the surface is 


& (testy) =O (6.3.3) 


hence 


Vas + Yus = Conant . (6.3.4) 


6.4 Modification of the Thermodynamic Fquation 


Writing (6.2.17) in terms of potential temperature gives: 


v falc) dO 
| Sibe; wig SPE 25 Sc ( 6.4.1) 


dee )andn( 3.72.3): 


(2. 
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Substitution of (6.4.2) into (6.4.1) and expansion of the total 


From equations 


derivative gives: 
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which is the thermodynamic equation modified for latent heat released 


during vertical motion. 


The following expression (which is a finite-difference 
approximation of the rate of change of the saturation mixing ratio 


for water vapour with temperature) is used to determine the effect 
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of the "correction" terms ( Gy and C uy ) on the thermodynamic equation: 


Ot AG 


where (Si = the change in the mixing ratio of liquid water resulting 


Ors; , Ate (6.4.4) 


from the application of the correction terms 


’ 


A® = temperature change at a given grid point resulting from 
the change in liquid water mixing ratio 
Yarg = saturation mixing ratio of water vapour 


Differentiating (6.2.9) and (6.2.10) gives: 


Stes z oe 7 (6.4.5) 
ky 


Since Nps ee equation (6.4.4) can be rewritten to obtain the 


expression for AS 


/ BG 
A8 = ieee ee 


‘ 
where A® = change in potential temperature per unit time 


After the inclusion of (6.4.6), equation (6.4.3) written in finite- 


difference form is: 
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The advection terms in this equation have been written for negative 


an n 
(Voting) and We. 


Changes in the equation of state and Poisson's equation 
resulting from the inclusion of water vapour and liquid water are 


ignored in this study. 
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CHAPTER '7 


THE DISSIPATION OF AN ALREADY EXISTING COASTAL FOG BANK BY 


THE SEA-BREEZE CIRCULATION 


ve inbroduction 

So far, only the results of a "dry" model have been presented. 
If the moisture equations of Chapter 6 are introduced into this study, 
the investigation of many "practical" meteorological situations becomes 
possible. This chapter shows the influence of the sea-breeze circulation 


on coastal fog. The fog is assumed to be already in existence at the 


initial time and its subsequent dissipation is studied in the model. 


Jee lnatiall Conditrone and Fields 


Thestactethat this model does notedesleexplicitly with solar 
radiation presents many problems which must be rectified by the use of 


rather drastic assumptions. 


The first problem that arises is the form of the temperature 
wave that is used to heat the land. Since the shape of the temperature 
wave depends on the amount of insolation reaching the ground, it is not 
possible to use Kuo's temperature wave because the fog reduces the 
amount of solar energy reaching the surface and also reduces the out- 
going long wave radiation. Hence, as a first approximation, a simple 


sine wave of the following form is used to heat the land: 


iver = ae + 5 Dm (15t) 


hcbuokacn, wind. site tebapemiaais 2a an an ) wae 
cue elit otal boculordat eo 3 awtieg e eebtenpe wrutehon ett IE 
panoserd ene ttmdia feo tyalenostas sao oe et Ye moaghtaw 
nolintier!n casorf=ase, ont to aohyulttet ade ewode weoqeds alt -eidteeog a 


a 


elt ts somiabee ih yoleald. 6 of Pedvegs et gol att «50? Lateugo 20. 
Lebon edt nt Selbute 2) softegieedh droupeaie ott bos =~ et 
* 9 - 
aloe tiv ylihe alge. Dest ster sb tiees ates sheds bi oa 
to ont ot ut eh ae sae nino nog . 
eS ae 


oa 
exons itt Sy eel Ag abe. Ast sateen ent at 
eurterequed wld To ayuda orth earth Boal, aft tod oF bogs et tadd evew 
ore ania mm me eae 


85 


erent ee = temperature of the land 
Vo 
e 


initial surface temperature 


i 


time in hours measured from 8:00 A.M. 


Geiger (1966) gives the variation of the sea temperature 
averaged over a depth of approximately one meter. The amplitude of the 
variation is about 0.1 C, is sinusoidal in character and reaches a 
maximum at 4:00 P.M. Since the amplitude above is the result of 
averaging over a depth of one meter, the temperature at the sea surface 
ig much larger. As a result, the following wave of amplitude 1.0 C is 


used to heat the sea: 


dcx 7, o.0G;9 * Av 115 (4-2) 


uhere Veeq = sea surface temperature 


tye = initial surface temperature 


+ = time in hours measured from 8:00 A.M. 


It will be assumed that none of the water vapour condensed 


falls but remains suspended in the air. 


The fog bank at the initial time extended from seaward lateral 
boundary to landward lateral boundary and was topped at a height of 550 
meters. The initial liquid water mixing ratio in the fog was set at 


1 gm/kg and zero elsewhere. 


The initial distribution of the water vapour mixing ratio was 


a function of height only: 
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eee) as era OS ZE 550m 
Anos 
re (zlos =e (550) > 550 £% £1950 m 


where Y,(2) = mixing ratio at height z 


Vy, = saturation mixing ratio 


This distribution gave saturation in the fog and a linear decrease of 


the mixing ratio above the fog. 


To comply with observations, the initial temperature distribution 
was chosen to give a weak inversion in the fog bank and a stable lapse 


rate above it: 


Bl et ©. Cole O£2£ 550 


T(z) = T(5s0) -.005 (2-550) 5504241950 
where \(z) = temperature at height z 


bls = initial surface temperature 


Using the same finite-difference parameters and number of grid 
points as in previous solutions, the initial fields of Table (7-1) were 


computed. 


7.3. Results for the Atmosphere Initially at Rest 


As a result of the assumptions mentioned above, it was felt 
that meaningful results would not be obtained after & hours. Results 
are given for the mid-point (12:00 noon) and the end of the period 


(4:00 P.M.). 
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The noon results (Figure (7-1)) showed that an offshore flow 
was present, resulting from the use of heat to evaporate the fog bank. 
The vertical velocity pattern showed weak subsidence over the sea and 
weak ascent over land. The fog bank over the sea was basically unchanged 
since the initial time while that over land had either cleared or "lifted" 


to, Lorm™ a stratus, cloud. 


By 4:00 P.M. (Figure (7-2)) a weak sea-breeze had formed and 
the vertical velocities had strengthened. The fog and stratus over land 
had completely dissipated. The fog was still present over the sea but 


its liquid water content was less than half its initial value. 


7.4 Results for an Initial Onshore Flow of 3 m/sec 


Lsmoelore, Only thcenoonwend 4: 00s PEM eresults ares presented, 


The noon results (Figure (7-3)) for this case showed that the 
fog had either cleared or lifted to stratus. Although the configuration 
of the fog bank over the sea was not significantly different from the 


initial time, its liquid water content had been greatly reduced. 


The 4:00 P.M. results (Figure (7-4)) showed that the fog had 
completely dissipated. Due to the increased forced convection resulting 
from the initial onshore wind, this case was more efficient in dissipating 


the fog bank than the case with no initial wind. 
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TABLE (7-1) Initial fields for fog dissipation study. 


av.6L5 vig OL ue Peg oi, « BR 2 
ee Se eee a 


1850 


1650 


1450 


1250 


1050 


850 


Height (M,) 


650 


450 


250 


50 


Height (M.) 


FIGURE (7-1) 


Fog 


—-40 


10 0 20 40 60 


Landward Distance from the Coast (Km.) 


FOG AND U COMPONENT 


——=— Fog -----U Component (m/sec) 


40 
(Km.) 
TEMPERATORE AND VemiivAe VERO y. 
--- Temperature (°C) 

— Vertical velocity (cm/sec) 


12:00 A.M. results for the atmosphere initially at 


20 
Lanaeyars Distance from the Coast 


&9 


TESts 


S).  .«- --- 4S eeleee 


90 


1850 


1650 


1450 


1250 


1050 


850 


Height (M.) 


650 


450 


250 


50 


-40 —20 0 20 40 60 
Landward Distance from the Coast (Km.) 


FOG AND U COMPONENT 


——-_ Fog -----U Component (m/sec) 


[OoUlt << -.s sae 


1650 
1450 
1250 
1050 


850 


Height (M.) 


650 


450 


250 


50 


60 


—40 —20 


20 4 
Landward Distance from the Coast (Km.) 


TEMPERATURE ANDO VER wICAL VELOC IY 
--- Temperature (°C) 
— Vertical velocity (cm/sec) 


FIGURE (7-2) 4:00 P.M. results for the atmosphere initially at rest. 
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FIGURE (7-3) 12:00 A.M. results for an onshore wind. 
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FIGURE (7-4) 4:00 P.M. results for an onshore wind. 
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CHAPTER 8 


CONCLUSIONS 


8.1 Contributions to Knowledge 


The original work contained in this study is the inclusion of 
the Coriolis term in the vertical equation of motion, the use of a non- 
hydrostatic mass conserving model to study synoptic control of the sea- 
breeze and the inclusion of equations for moisture to study a practical 


meteorological problem. 


8.2 Validity of the Model 


In spite of all the assumptions made, a comparison of the 
results of this study with observations shows that the sea-breeze is 
well reproduced by this model. The greatest difference between the 
theoretical results and the observations is in the vertical extent of 
the sea-breeze. Koschmieder and Hornikel (1942) observed that the 
vertical penetration of the sea-breeze under calm gradiemt wind 
conditions was less than that under offshore gradient wind conditions. 
However, this study predicts the opposite. Table (8-1) summarizes the 


model's performance in representing the available observations. 


§.3 Suggestions for Future Study 


The energy equation should be used in each of the sea-breeze 


cases studied in order to fully understand the circulation. 


i ee 


i: a 
as 


wis 


to nokautait ood ot Giada and) At hemtegdon | 
~soe 5 to sew edt notion Yo apltaupe Teoksrev edt at 
~ese oii ‘to foatsies obtyoree Yitnhe ot lekee | 


| aout 
one i a 


i ee of A nly 7 
as a a a ee 
+ = ties abcabenanae Weome mS, 
i susond—nae aid dats awciin entpriindé date ybatre why Yo 29 ie 
wid useoted sanmislteb Jewtaety ent cet se aaa 
nis duibixs Sentvnal add ot pt anectiaaedo ant hea qhtinen 
edit fait terigedo (SRCE) Dettmndll baw iebatmriseat mori 
fiw dentin aso dobel MaNBA-aes oir Io aoctendeeey Bp 
-anststinos Bie Sanchar StonitTe nob Sar amy wnat 2 
eer ee, RAE 


oy 


calm gradient 
wind 


offshore gradient 
wind case 


onshore gradient 
wind case 


fog dissipation 


OBSERVED FEATURES 
DESCRIBED BY MODEL 


hodograph 


spreading out of 
wind from coast— 
line 


weak reverse flow 
allot t 


beginning at sea 


front. well 
Lepresented 


less landward 
penetration than 
calm wind case 


sea-breeze weak 


average increase 
OL twain m/sec, 
model gives 1.5 
m/sec 


no observations 
available 


TABLE (8-1) Performance of the model. 
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In the cases where a synoptic wind is used, the model should 
be allowed to reach a steady state, i.e. the synoptic wind should be 
allowed to produce an Ekman spiral before the temperature wave is 


applied to the land. 


The most significant improvement that could be made is the 


inclusion of radiation in the model. 


The model could also be improved by more accurate equations 
for the lower sub-layer. More realistic results could be obtained if 
Chorin's method is modified to use a compressible equation of continuity. 
With the shallow convection restriction removed, the number of 
meteorological problems that could be investigated is greatly increased, 


such as the enhancement of cumulus development by the sea-breeze. 
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APPENDIX A 


The following is a listing of the computer program used 
iNeEneet oe dissipation study of Chapter 7. 


SUBROUTINE MAPPA is courtesy of Y. Mahrer. 
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